RNA polymerase motors: dwell time distribution, velocity and dynamical phases 
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Polymerization of RNA from a template DNA is carried out by a molecular machine called RNA 
polymerase (RNAP). It also uses the template as a track on which it moves as a motor utilizing 
chemical energy input. The time it spends at each successive monomer of DNA is random; we 
derive the exact distribution of these "dwell times" in our model. The inverse of the mean dwell time 
satisfies a Michaelis-Menten-like equation and is also consistent with a general formula derived earlier 
by Fisher and Kolomeisky for molecular motors with unbranched mechano- chemical cycles. Often 
many RNAP motors move simultaneously on the same track. Incorporating the steric interactions 
among the RNAPs in our model, we also plot the three-dimensional phase diagram of our model for 
RNAP traffic using an extremum current hypothesis. 

PACS numbers: 



I. INTRODUCTION 

RNA polymerase (RNAP) is a molecular motor \[^. It 
moves on a stretch of DNA, utilizing chemical energy in- 
put, while polymerizing a messenger RNA (mRNA) [2]. 
The sequence of monomeric subunits of the mRNA is 
dictated by the corresponding sequence on the template 
DNA. This process of template-dictated polymerization 
of RNA is usually referred to as transcription [3] . It com- 
prises three stages, namely, initiation, elongation of the 
mRNA and termination. 

We first report analytical results on the characteris- 
tic properties of single RNAP motors. In our approach 
P], each RNAP is represented by a hard rod while the 
DNA track is modelled as a one-dimensional lattice whose 
sites represent a nucleotide, the monomeric subunits of 
the DNA. The mechano-chemistry of individual RNAP 
motors is captured in this model by assigning m distinct 
"chemical" states to each RNAP and postulating the na- 
ture of the transitions between these states. The dwell 
time of an RNAP at successive monomers of the DNA 
template is a random variable; its distribution charac- 
terizes the stochastic nature of the movement of RNAP 
motors. We derive the exact analytical expression for the 
dwell-time distribution of the RNAPs in this model. 

We also report results on the collective movements of 
the RNAPs. Often many RNAPs move simultaneously 
on the same DNA track; because of superficial similari- 
ties with vehicular traffic [5,], we refer to such collective 
movements of RNAPs as RNAP traffic 0, 0, 0, S, Q- 
Our model of RNAP traffic can be regarded as an exten- 
sion of the totally asymmetric simple exclusion process 
(TASEP) [10] for hard rods where each rod can exist at 
a location in one of its m possible chemical states. The 
movement of an RNAP on its DNA track is coupled to 
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the elongation of the mRNA chain that it synthesizes. 
Naturally, the rate of its forward movement depends on 
the availability of the monomeric subunits of the mRNA 
and the associated "chemical" transitions on the domi- 
nant pathway in its mechano-chemical cycle. Because of 
the incorporation of the mechano-chemical cycles of in- 
dividual RNAP motors, the number of rate constants in 
this model is higher than that in a TASEP for hard rods. 
Consequently, we plot the phase diagrams of our model 
not in a two-dimensionl plane (as is customary for the 
TASEP), but in a 3-dimensional space where the addi- 
tional dimension corresponds to the concentration of the 
monomeric subunits of the mRNA. 



II. MODEL 

We take the DNA template as a one dimensional lat- 
tice of length L and each RNAP is taken as a hard rod of 
length I in units of the length of a nucleotide. Although 
an RNAP covers i nucleotides, its position is denoted 
by the leftmost nucleotide covered by it. Transcription 
initiation and termination steps are taken into account 
by the rate constants uOqc and a;/3, respectively. A hard 
rod, representing an mRNA, attaches to the first site 
i = 1 on the lattice with rate uj^ if the first I sites are 
not covered by any other RNAP at that instant of time. 
Similarly, an mRNA bound to the rightmost site i = L is 
released from the system, with rate uof^. We have assumed 
hard core steric interaction among the RNAPs; therefore, 
no site can be simultaneously covered by more than one 
RNAP. At every lattice site i, an RNAP can exist in one 
of two possible chemical states: in one of these it is bound 
with a pyrophosphate (which is one of the byproducts of 
RNA elongation reaction and is denoted by the symbol 
PPi), whereas no PPi is bound to it in the other chemical 
state (see figH]). For plotting our results, we have used 
throughout this paper 0012 = 31.4 s~^,iJi2 — 30.0 
and uj 21 = [NTP] Co 21 where [NTP] is concentra- 
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tion of nucleotide triphosphate monomers (fuel for tran- 
scription elongation) and 0J21 = 10^ M~^s~^. 
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FIG. 1: The mechano- chemical cycle of a RNAP in our model. 
In the chemical state 1 the RNAP is bound to pyrophosphate 
(PPi) whereas no PPi is bound to it in the state 2. The 
integer subscript j labels the position of the RNAP motor on 
its track. 



III. DWELL TIME DISTRIBUTION 

For every RNAP, the dwell time is measured by an 
imaginary "stop watch" which is reset to zero whenever 
the RNAP reaches the chemical state 2, for the first time^ 
after arriving at a new site (say, i + 1-th from the i-th). 

Let be the probability of finding a RNAP in the 
chemical state /i at time t. The time evolution of the 
probabilities P^ are given by 
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There is a close formal similarity between the mechano- 
chemical cycle of an RNAP in our model (see figH]) and 
the catalytic cycle of an enzyme in the Michaelis-Menten 
scenario The states 2 and 1 in the former corre- 

spond to the states E and ES in the latter where E rep- 
resents the free enzyme while ES represents the enzyme- 
substrate complex. Following the steps of calculation 
used earlier by Kuo et al. [12] for the kinetics of single- 
molecule enzymatic reactions, we obtain the dwell time 
distribution 



fit) = i^Mi^ - m - exp{-{A ■ 



where 



2A 



A 



B)t}] (3) 



(^12 ■ 



'12 ■ 



/ j2 



^21 



■ ^12^21 



B 



L0\2 + <^12 ■ 



21 



(4) 



(5) 



O 

qI 



0.7 
0.6 \ 
0.5 
0.4 
0.3 
0.2 
0.1 




[NTP]=0.5 \M 
[NTP]=1.0 ]M 



CO12 = 20 s 



2 4 6 

Dwell time (s) 



10 12 



(b) 



O 

qI 



0.7 
0.6 \ 
0.5 
0.4 
0.3 
0.2 
0.1 




CO1 2=5.0 s"^ 
CO1 2=20.0 s'"* 



[NTP] = 1 .0 laM 



4 6 
Dwell time (s) 



10 12 



FIG. 2: The dwell time distribution of a single RNAP is our 
model is plotted for (a) two different values of NTP concentra- 
tion, keeping the rate of PP^ release (and, hence, a;i2) fixed, 
and (b) two different values of a;i2, keeping the NTP concen- 
tration fixed. 



The dwell time distribution (|3|) is plotted in figiS De- 
pending on the magnitudes of the rate constants the peak 
of the distribution may appear at such a small t that it 
may not possible to detect the existence of this maxium 
in a laboratory experiment. In that case, the dwell time 
distribution would appear to be purely a single expo- 
nential [1^, [3, [iBj. It is worth pointing out that our 
model does not incorporate backtracking of RNAP mo- 
tors which have been observed in the in-vitro experiments 
[H, • It has been argued by some groups [18] that 
short transcriptional pausing is distinct from the long 
pauses which arise from backtracking. In contrast, some 
other groups [19] claim that polymerase backtracking can 
account for both the short and long pauses. Thus, the 
the role of backtracking in the pause distribution remain 
controversial. Moreover, it has been demonstrated that 
a polymerase stalled by backtracking can be re-activated 
by the "push" of another closely following it from behind 
[20]. Therefore, in the crowded molecular environment of 
intracellular space, the occurrence of backtracking may 
be far less frequent that those observed under in-vitro 
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conditions. Our model, which does not ahow backtrack- 
ing, predicts a dweh time distribution which is quahta- 
tively very similar to that of the short pauses provided 
the most probable dwell time is shorter than 1 s. 

From equation ([3|) we get the inverse mean dwell time 
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where Vmax = ^12 and Km = (^12 + ^i2)/^2i- The 
form of the expression (|6]) is identical to the Michaelis- 
Menten formula for the average rate of an enzymatic 
reaction. It describes the slowing down of the "bare" 
elongation progress of an RNAP due to the NTP reac- 
tion cycle that it has to undergo. The unit of velocity is 
nucleotide / second. 

The fluctuations of the dwell time can be computed 
from the second moment 
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of the dwell time distribution. We find the randomness 
parameter [H, |2l|, [22| 



by dlO]), is plotted against the NTP concentration in 
fig El for three different values of CJ12. At sufficiently low 
NTP concentration, 7 is unity because NTP binding with 
the RNAP is the rate-limiting step. As NTP concentra- 
tion increases, 7 exhibits a nonmonotonic variation. At 
sufficiently high NTP concentration, PP^-release (which 
occurs with the rate UJ12) is the rate-limiting step and, 
therefore, 7 is unity also in this limit. This interpretation 
is consistent with the fact that the smaller is the mag- 
nitude of a;i2, the quicker is the crossover to the value 
7 = 1 as the NTP concentration is increased. 

The randomness parameter yields the diffusion coeffi- 
cient [21] 
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The expression (p!2]) is in agreement with the general ex- 
pression for the effective diffusion constant of a molecular 
motor with unbranched mechano-chemical cycle which 
was first reported by Fisher and Kolomeisky [23]. 



IV. PHASE DIAGRAMS 
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Now we will take into account the hard core steric in- 
teraction among the RNAPs which are simultaneously 
moving on the same DNA track. Equations ([T] and [2]) 
will be modified to 
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FIG. 3: The randomness parameter 7 is plotted against NTP 
concentration for three values of the parameter 0012. 

Note that, for a one-step Poisson process f{t) = 
ix}exp{—ix}t)^ 7 = 1. The randomness parameter 7, given 



= 0012 Pi - 4l ^2 g + C^i'2 Pi Q (14) 

Where Q is conditional probability [?] of finding site 

{i — i — 1 for backward motion) vacant, given there is a 

particle at site i. 

Due to the steric interactions between RNAP's their 
stationary flux J (and hence the transcription rate) is 
no longer limited solely by the initiation and release at 
the terminal sites of the template DNA. We calculate the 
resulting phase diagram utilizing the extremum current 
hypothesis (ECH) 0, |25|. The ECH relates the flux 
in the system under open boundary conditions (OBC) to 
that under periodic boundary conditions (PBC) with the 
same bulk dynamics. In this approach, one imagines that 
initiation and termination sites are connected to two sep- 
arate reservoirs where the number densities of particles 
are p- and p+ respectively, and where the particles fol- 
low the same dynamics as in the bulk of the real physical 
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system. Then 



flux is given by the expression 
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The actual rates cOa and cof^ of initiation and termination 
of mRNA polymerization are incorporated by appropri- 
ate choice of /)_ and respectively. 
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FIG. 4: The 3-d phase diagram of our model for RNAP traffic. 
The LD and HD phases coexist on the surface /. The surfaces 
// and /// separate the MC phase from the HD and LD 
phases, respectively. 




FIG. 5: Projections of several 2-d phase diagram of RNAP 
traffic on the iJa-(^f3 plane for several values of i02i- The num- 
bers on the phase boundary lines represent the value of 
The inclined lines have LD and HD above and below, respec- 
tively, while the MC phase lies in the upper right corner. 



By comparing (p!6|) with the exact current-density rela- 
tion of the usual TASEP for extended particles of size 
i [2^, [23, d^, which have no internal states (formally 
obtained by taking the limit UJ12 ^ 00 in the present 
model), we predict that the stationary current (i.e. the 
collective average rate of translation) is reduced by the 
occurrence of the intermediate state 1 through which the 
RNAPs have to pass. 

From (p!5|) one expects three phases, viz. a maximal- 
current(MC) phase with with bulk density p*, a low- 
density phase (LD) with bulk density p_, and a high- 
density phase (HD) with bulk density Using argu- 
ments similar to those used in Ref. [29] in a similar con- 
text, we get ^] 
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An expression for J{p) was reported by us in Ref. Q. 
In the special case where the dominant pathway is that 
shown in Fig. [T] {uji2 = foi" further calculation as 
U12 ^ ^12)5 have 
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FIG. 6: Same as in fig|5] except that the projections are on 
(jtj 21-00 f3 plane for several values of Ua. The inclined lines have 
LD and HD above and below, respectively. Each vertical line 
separates the LD phase on the left from the MC phase on its 
right. 

The condition for the coexistence of the high density 
(HD) and low density (LD) phases is 



J (p.) = j{p+) 



(20) 
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with p_ < p* < p+. Using the expression ([T6|) for J in 
(f2Q|) we get 
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Substituting (p!8|) and (p!9|) into ([2T]) , we get the equation 
for the plane of coexistence of LD and HD to be uof^ = 
f{ua,uj^i) where 
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In order to compare our result with the 2-d phase dia- 
gram of the TASEP in the — C(;/3-plane, we project 2-d 
cross sections of the 3-d phase diagram, for several dif- 
ferent values of onto the cja — oop-plane. The lines of 
coexistence of the LD and HD phases on this projected 
two-dimensional plane are curved , a similar curvature 
is also reported by Antal and Schiitz [30]. This is in 
contrast to the straight coexistence line for LD and HD 
phases of TASEP. 



The bulk density of the system is guided by following 
{^^ations: 
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In Fig. m we plot the 3d phase diagram. 




FIG. 7: Same as in fig|5] except that the projections are on 
00^1-00 a plane for several values of cj/^. Here the inclined lines 
have HD and LD, respectively, above and below. Each vertical 
line separates the HD phase on the left from the MC phase 
on its right. 

In general, a plane 6^21= const ant intersects the surfaces 
I, II and HI thereby generating the phase transition lines 
between the LD, HD and MC phases in the 0Oa-^f3 plane. 
We have projected several of these 2d phase diagrams, 
each for one constant value of LO21 in figure \5\ In the 
inset, we have shown the value of UJ21 for different lines. 



We have also projected several 2d phase diagrams in the 
0021-00(3 plane and 00 21-00^3 plane, respectively, in figures [6] 
and [71 



V. SUMMARY AND CONCLUSION 

In this paper we have reported the exact dwell time 
distribution for a simple 2-state model of RNAP motors. 
From this distribution we have also computed the average 
velocity and the fluctuations of position and dwell time 
of RNAP's on the DNA nucleotides. These expressions 
are consistent with a general formula derived earlier by 
Fisher and Kolomeisky for a generic model of molecular 
motors with unbranched mechano-chemical cycles. 

Taking into account the presence of steric interactions 
between different RNAP moving along the same DNA 
template we have plotted the full 3d phase diagram of 
a model for multiple RNAP traffic. This model is a bi- 
ologically motivated extension of the TASEP, the novel 
feature being the incorporation of the mechano-chemical 
cycle of the RNAP into the dynamics of the transcrip- 
tion process. This leads to a hopping process with a dwell 
time distribution that is not a simple exponential. Nev- 
ertheless, the phase diagram is demonstrated to follow 
the extremal-current hypothesis [IJ] for driven diffusive 
systems. Using mean field theory we have computed the 
effective boundary densities that enter the ECH from the 
reaction constants of our model. We observe that the 
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collective average rate of translation as given by the sta- 
tionary RNAP current ([16]) is reduced by the need of the 
RNAP to go through the pyrophosphate bound state. 
This is a prediction that is open to experimental test. 

The 2d cross sections of this phase diagram have 
been compared and contrasted with the phase diagram 
for the TASEP. Unlike in the TASEP, the coexistence 
line between low- and high-density phase is curved for 
all parameter values. This is a signature of broken 
particle- vacancy symmetry of the RNAP dynamics. The 
presence of this coexistence line suggests the occurrence 



of RNAP "traffic jams" that our model predicts to 
appear when stationary initiation and release of RNAP 
at the terminal sites of the DNA track are able to 
balance each other. This traffic jam would perform an 
unbiased random motion, as argued earlier on general 
theoretical grounds in the context of protein synthesis 
by ribosomes from mRNA templates |31l] . 
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